Two parameter flow of a xx {u) — a xy (u) for the graphene quantum Hall system in ac 

regime 
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Flow diagram of (a xx ,a xy ) in finite-frequency (u) regime is numerically studied for graphene 
quantum Hall effect (QHE) system. The ac flow diagrams turn out to show qualitatively similar be- 
havior as the dc flow diagrams, which can be understood that the dynamical length scale determined 
by the frequency poses a relevant cutoff for the renormalization flow. Then the two parameter flow 
is discussed in terms of the dynamical scaling theory. We also discuss the larger-o; regime which 
exhibits classical flows driven by the raw frequency uj. 



I. INTRODUCTION 

In the quantum Hall effect (QHE), one standard and 
graphically clear way to grasp the physics involving the 
localization effect is the a xx — a xy diagram, in which 
we look at the scaling flow (trajectories when the sam- 
ple size is varied) of the longitudinal conductivity and 
Hall conductivity (<7 XX , a xy ). The scaling property of 
the static QHE system, especially the quantization of 
the Hall conductivity into the multiples of e 2 /h and van- 
ishing longitudinal conductivity, is beautifully captured 
with the a xx — a xy diagram as originally discussed by 
Pruisken and Khmclnitskii in terms of the non-linear 
sigma modeli^— . For the conventional two-dimensional 
electron gas (2DEG), there exist (i) stable fixed points at 
(p xx ,a xy ) = (0,n)(n: integer), along with (ii) unstable 
fixed points characterizing derealization at {<r xx ,cr xy ) = 
{<J XX , n + 1/2). The attraction into the former, quantum- 
Hall fixed point accounts for the Hall insulating state 
with quantized values of Hall conductivity, while the lat- 
ter, unstable fixed point accounts for delocalized states 
at the center of each Landau level (LL), and dominates 
the behavior of the plateau-to-plateau transition in a xy . 

Scaling properties of the Anderson transition have at- 
tracted both theoretical and experimental interests, since 
they should be universal and depend only on the symme- 
try class of the system^. The critical exponent has been 
numerically studied for the lowest LL£, and later with the 
Chalkcr-Coddington network model 6 - and experimentally 
confirmed by Li et al^ The universal value of longitudinal 
conductances at the LL centers are intensively discussed 
with a tight binding lattice model^ . Thus the scaling 
behavior at the plateau-to-platcau transition has been 
established. 

On the other hand, rapid advances in the terahertz 
(THz) spectroscopy technique have made the optical 
responses of the quantum Hall system, such as cy- 
clotron resonances and Faraday rotations, experimentally 
accessibl e 9 ' 10 . Specifically, the Faraday rotation is pro- 
portional to the optical Hall conductivity a xy (uj), and we 
have an intriguing problem of how the static Hall conduc- 
tivity, which may be regarded as a topological quantity-ii, 
evolves into the optical Hall conductivity, especially in 
the relevant (cyclotron) energy scale which falls upon 



the THz regime*^— On the optical Hall conductivity, 
we have recently shown, numerically, that the plateau 
structure in a xy (uj) is unexpectedly retained in the ac 
(THz) regime in both 2DEG and in graphene, although 
the plateau height deviates from the quantized values in 
acJ^. Graphene is particularly interesting, since a mass- 
less Dirac system is realized as the low-energy physics, 
which accommodates a novel Dirac QHE is observe d 16 ' 17 , 
for which the scaling theory of QHE in graphene has been 
formulated in terms of the non-linear sigma model^. Ex- 
perimentally, the ac plateau has been observed in a GHz 
Faraday rotation measurement for a 2DEG system-^ and 
recently in THz regime 2 ^. In the graphene QHE system, 
we expect plateau structures at the tail (small frequency) 
region, while Faraday rotation was measured in the re- 
gion around cyclotron resonances 2 ^. 

Considering these advances in spectroscopies of QHE 
systems, it is important to study the systematic behav- 
ior of plateau structure in ac and optical regimes. Ro- 
bustness of the ac plateau structure against disorder as 
revealed in the numerical result can be understood if we 
consider the effect of localization, which dominates the 
physics of electrons around the centers of Landau levels 
in disordered QHE systems, in finite-frequency regime, 
following the scaling theory of the Anderson transition. 
Namely, a finite frequency put an effective cutoff (the dy- 
namical length scale L u ) for the system, and the plateau 
in the ac Hall conductivity should be retained in the re- 
gion where the localization length (£ ~ \e — £ c \~ u ) di- 
verging toward LL center is smaller than the dynamically 
posed cutoff L u . In the two parameter flow picture, this 
can be viewed as the dynamical length scale L u deter- 
mining a scale where the renormalization is stopped. 

The dynamical scaling behavior has been studied for 
the ac longitudinal conductivity^, and for the optical 
Hall a xy (uj)Q. Now, it is interesting to combine both 
the optical longitudinal and Hall a xy (iv) and numerically 
map out the two parameter (a xx — o xy ) flow diagram in 
the ac regime. 

In the present work we have calculated both optical 
longitudinal and Hall conductivities (a xx , a xy ) for the 
graphene QHE system with a potential disorder, and 
combined them to numerically examine the two parame- 
ter a xx — <j xy flow diagram in the ac regime. We study 
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n = LL in the graphcnc QHE system with exact diag- 
onalization method to treat disorder effects. One partic- 
ular point of interest is the behavior around fixed points 
in the a xx (uj) — a xy (uj) diagram. There, we have focused 
on the n = Dirac Landau level, where the peculiar- 
ity of graphene appears as the property that n = is 
an electron-hole symmetric point. In a small-cj regime, 
we obtain numerical results which are coherent with the 
above picture that the <j xx (uj) — a xy (uj) flow obeys the 
Pruisken's two-parameter flow, with as a relevant 
cutoff for the system in ac region, where the flows are 
between a xy = ±2e 2 /h reflecting the graphene QHE in- 
cluding the valley and the spin degeneracies. We also dis- 
cuss a large-w regime where the frequency uj is compara- 
ble with the cyclotron frequency uj c and exhibits classical 
flows driven by the raw frequency 



related to L with = L 2 /2ir£ 2 in a finite system. From 
the eigenfunctions ip a and eigenenergies e a obtained with 
the exact diagonalization, the optical Hall conductivity^ 
is evaluated from the Kubo formula 2 ^ as 
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where /(e) is the Fermi distribution, and r\ a low-energy 
cutoff. The current matrix element, j x b , has a selection 
rule peculiar to Dirac model (n O ±n ± 1 with n the 
Landau index), which is distinct from that (no n ± 1) 
for 2DEG as 

3x' n = ev F C n C n > [sgn(n)S M _ 1 ^ n ,\ + sgn(n')<5 H+lj |„,|] , 
jy' n =iev F C n C n < [sgn(n)6\ n \_ lt \ n ,\ - sgn(7i')£|n|+i,|rt'|] , 



II. FORMALISM 

For graphene QHE system, we employ the two- 
dimensional effective Dirac model, 



H 



vpT z er ■ n 



V(r), 



(1) 



where vp is Fermi velocity, <r = (a x , a y ) and t z the Pauli 
matrices acting on the space of two sublattices (A, B) 
and two valleys (K, K'), tv = p + eA with p = (p x ,p y ) 
the momentum, and A the vector potential. Disorder is 
introduced by a random potential, 

V{r) = exp(-|r - R u \ 2 /2d 2 )/ (2ird 2 ), 



composed of Gaussian scattering centers of range d and 
Ui.j takes a value in (— u,u) randomly. Here we take 
d = 0.7£, where £ = yJh/eB is the magnetic length. 
For numerical facility, impurity sites Rij are periodically 
placed on Rij = (2ir£ 2 / L)(i, j) with L being the linear 
dimension of the sample. A measure of disorder is given 
by the Landau level broadening^, T = 2u[Ni mp /2n(£ 2 + 
2d 2 )L 2 ] 1 / 2 , with a number of impurity sites Ni mp . We 
assume smooth potential disorders in the length scale of 
underlying lattice structure, and we neglect inter-valley 
scattering. The cyclotron energy is, for a Dirac particle, 
given by lu c = \/2vf/£- 

Since we are interested in the dynamical <j xx (uj) — 
<j xy (uj), which should be related to the localization 
physics, we obtain the eigenstates of the Hamiltonian 
with an exact diagonalization, which is done for a sub- 
space spanned by a finite number of Landau levels (LL's) 
around n = LL, for L x L systems with L/£ varied over 
20,30,40. Here we retain 5 LLs (n = -2 ~ 2), which 
poses an ultraviolet cutoff.— In the Landau gauge A ~ 
(0,Bx), the basis function is ip n ,k = e~ lky (j) n (x — £ 2 k y ), 
where 4> n is the Dirac-Landau function in the n-th Lan- 
dau level^, and wavcnumbers k takes an integer mul- 
tiples of 2tt/L with a periodic boundary condition for 
2/-dircction. The number of discrete wavcnumbers Nk is 



where C„ = 1 (n = 0) or 1/ \[2 (otherwise^i^l. 

The longitudinal conductivity, on the other hand, is 
given by 

_ h ^ f(e b ) - f(e a ) \jf\ 2 V 
xx (UJ) ~ L 2 £b - £a ( E5 _ Ea _ M 2 + V 2 ■ 

We note that the low-energy cutoff r), which affects the 
uj ~ behavior of <j xx (uj), should be chosen close to the 
Thoulcss energy, which is typically of the order of the 
energy level spacing ~ l/L 2 ^— The temperature in 
the Fermi distribution function /(e) is set to be small as 
far as the low- frequency behavior of a xx (oj) is numerically 
stable, which is achieved if we put the temperature to be 
of the order of the level spacing ~ 1/L 2 . For the scaling 
analysis the calculation is repeated for varied sample size 
L, Fermi energy £p and frequency uj. Throughout the 
paper the length, energy and frequency are respectively 
in units of £, huj r and uj r . 



III. 



DIAGRAM IN AC REGIME 



First we discuss the behavior a xx (oj) and a xy (uj) sep- 
arately. <j xx (eF,w) in FigQJa) shows a ridge struc- 
ture along the uj axis when the Fermi energy tp is 
around the delocalized region at each LL. When we 
increase L in FigfTJb), the width of a xx (oj = 0) be- 
comes narrower, while the peak height stays almost con- 
stant which is expected from the universal longitudi- 
nal conductance^ ' 31 ' 32 . On the other hand, a xy (eF,u) 
plotted against ep in Figfljc) shows a transition from 
a xy = — 2 plateau to a xy — 2 plateau around n = LL 
with 2 valley and 2 spin degeneracies, where the transi- 
tion width of <i xy {uj = 0) sharpens with increasing L. 

Now we are in a position to examine a xx — a xy dia- 
gram in Figj2] First, the diagram for uj = is depicted 
in FiglUJa), where almost all the points arc attracted 
to the points (a xx , o xy ) = (0, ±2)e 2 //i with increased 
sample size L, while the point sitting at a xy = only 
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FIG. 1: (a) <j xx (eF,oj) plotted against the Fermi energy 
€f and the frequency u for L = 30. Lower panels depict 
o~ xx (u) = 0) (b) and <J xy (u> = 0) (c) for various sample sizes 
with disorder strength T /Tiu c = 0.4. 



FIG. 2: The flow of (<j xx (lj) — a xy (u})) in the graphene quan- 
tum Hall system with a disorder strength T/hcj c = 0.4. (a) 
The flow in dc regime (uj — 0) for various Fermi energy ef 
and system size L/i = 20,30,40. (b) The flow in ac regime 
for various values of Fermi energy e f and frequency uj with a 
fixed system size Ljt — 30 



exhibits a tiny, upward flow. This numerically calcu- 
lated dc flow diagram is clearly understood in terms of 
Pruisken's two parameter flow picture^— . We can inter- 
pret the static result that the flows starting from a xy ^ 
corresponding to those flowing into the stable fixed points 
at (<r xx ,a xy ) = (0,±2)e 2 //i which describes Hall plateau 
for graphene, while the tiny upward flow around a xy = 
corresponds to the flow that starts from the point with 
a <r xx smaller than a xx (rather than larger a xx expected 
from SCBA values), and shows a renormalization to the 
unstable fixed point at (a xx ,a xy ) = (a xx ,0). So we are 
seeing the region below the unstable fixed point. This 
upward flow toward the unstable fixed point would re- 
flect the existence of the delocalized state at the LL cen- 
ter, since it is percolating through the sample and has 
a metallic nature. The longitudinal conductivity then 
increases with the sample size and converges to the uni- 
versal conductance at the LL center. This dc result for 
graphene n — LL is consistent with Nomura ct al who 
discusses the Thouless-number and the Hall conductivity 
for the dc flow diagram^ 

If we now turn to the ac result that is plotted for a fixed 
system size L/£ = 30 varying frequency uj /uj c = 0.0025 ~ 
0.015 in Figl^b), where a behavior quite similar to the dc 



data is found. Namely, almost all the points away from 
the LL center are attracted to the QHE fixed points, 
while the point on a xy = at the center of LL flows 
only slightly shift upwards. This behavior is understood 
that, in this small frequency regime, the relevant cutoff 
length scale for the critical behavior of localization length 
£ or the renormalization equation for two-parameter flow 
is posed by the frequency through the dynamical length 
scale L^j instead of the sample size L in the dc regime, 
so that the overall behavior is determined by the same 
two-parameter flow, where the effective cutoff alone is 
changed systematically with the frequency u as L u ~ 



IV. DYNAMICAL SCALING ANALYSIS 

Now let us describe the dynamical scaling for 
&xx{£f,w) and a xy {eF,uS). The scaling argument starts 
from an ansatz that the optical conductivity depends on 
Fermi energy sp and frequency uj only through the ra- 
tios L/£ and L u /^. The physical quantities should then 
be described in terms of the universal scaling function of 
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FIG. 3: (a) The flow of (a xx (u>) — <j xy (u)) in the graphene 
quantum Hall system for renormalized frequency ujL z = 3 
with z = 2 and L = 20, 30, 40 with a disorder strength 
T/hio c — 0.4. (b) Flows when the sample size is varied as 
L = 20 — > 40 for a fixed of toL z . For each value of uL z we 
plot the flows corresponding to various values of Fermi energy 
e F . The value of is indicated by different symbols (circles, 
squares, etc) that mark the smallest sample size (L = 20). 
The results for various values (color-coded) of ujL z =3 — 6 
are superposed. 



the ratios L/£ and Here £ is the localization length 

with a critical behavior £ ~ l/\s F — e c \ v , where e c is the 
critical energy which coincides with the center of the LL 
(e c = for n = 0), and v the localization critical expo- 
nent. The dynamical length scale, which is the distance 
over which an electron travels during one cycle, l/u>, of 
the ac field, is assumed to behave as L u ~ l/w 1 / 2 , where 
z is the dynamical critical exponent, assuming z = 2 in 
this paper since we treat the non-interacting electrons^, 
while z — 1 is established for the case with electron- 
electron interaction^. For these critical behaviors the 
dynamical scaling ansatz for the longitudinal and trans- 
verse conductivities amounts to2S 



a xx (£ F ,u),L) = —F xx {Se F L 1/v ,ujL z ), 
e 2 

o xy (e F , u, L) = —FxyiSepL 1 ^ \ojL z ), (3) 



where F xx , F xy are universal scaling functions, and 5e F = 
e F — e c . 

This ansatz is only valid for the critical region, where 
the deviation of Fermi energy from the LL center 5e F 
is assumed to be small and the frequency oj also to be 
small. In Sec. 3 and Sec. 4, we precisely consider this re- 
gion with small Se F and oj, while in Sec. 5 where we shall 
treat a large-frequency region this ansatz should be no 
longer applicable. The ansatz indicates that the flow 
of (<j xx (ijj),<j X y(u))) in the ac region depends on the fre- 
quency oj only through the ratio of the rescaled frequency 
with the system size as u>L z . 

We now interpret the behavior of a xx — a xy diagram 
in terms of the dynamical scaling eqn.Q as shown in 
Figj3] When we consider a dynamical scaling behavior, 
it is convenient to consider a rescaled frequency ujL z since 
a dependence on the frequency appears in a form of ojL z 
in eqn.([3|). In Fig. [3ja) wc show a result for a fixed 
rescaled frequency loL z = 3, then the flow of ac con- 
ductivities can be discussed in terms of varying L as in 
dc case. The dynamical scaling hypothesis (eqn.Q) ex- 
pects that (a xx ,a xy ) right on e F = with itpL 1 ^ = 
for all L should depend only on wL z , i.e., does not flow, 
while at the center of Figj3ja) a slight upward flow is seen 
showing a metallic behavior discussed above. Away from 
e F = 0, on the other hand, the flow should depend only 
on SepL 1 ^ from eq.([3]) for a fixed value of the rescaled 
frequency loL z , which implies that the flows starting from 
various e F ^ should reside on a single curve as seen in 
FigSta). 

In order to examine the dependence of the two param- 
eter flow on the rescaled frequency ojL z , we superpose 
all the results in Fig[3]Jb). There, we show flows when 
the sample size is varied as L = 20 — > 40 for fixed values 
of ojL z . For each value of u>L z we plot the flows cor- 
responding to various values of Fermi energy e F . The 
value of e F is indicated by different symbols that mark 
the smallest sample size (L = 20). Fig[3]Jb) thus visual- 
izes the the two-parameter flows superposes for various 
values of cuL 2 and for various values of Fermi energy e F . 
The results for various values of ujL z = 3 — 6 are then 
superposed in FigJHlJb). In this summary plot we can see 
that the (a xx , a xy ) flows for different values of u>L z tend 
to coalesce into a single curve in the region away from 
&xy = 0. This is a consequence that we see the same two 
parameter flow with various cutoff length scale posed by 
different rescaled frequencies. Close to the unstable fixed 
point at {(J xx ,a xy ) = (<J XX ,0) the flow shows a metallic 
behavior (i.e., a xx increasing with L) renormalizing into 
the unstable fixed point with increasing a xx , and slightly 
deviates from a single curve. The role of the value of 
increasing ojL z appears in a shift of the initial position of 
the flow in (cr xx (uj), a x y(uj)) for each e F toward the op- 
posite direction to the flow. More precisely, for a larger 
ujL z we have more broadened peak structure in a xx and 
more broadened transition width in a xyi so that the ini- 
tial point (smallest-L data) of the flow for each e F shifts 
closer to the unstable fixed point (the center of the flow 
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FIG. 4: (a) cr X x(eF = 0, uj) plotted for frequency up to uj — 
0.1, where the data for various system sizes almost coincide 



with each other for large u). (b) a xx {ui) 
raw frequency oj with L = 30. 



a X y{ui) diagram for 



at a xy ) = (o"£ x , 0) ), which we can observe in a shift 
of initial values in FigJ3Jb). This behavior arises from the 
fact that the relevant cutoff for the two-parameter flow 
in the ac regime is determined by the dynamical length 
scale L u , where a larger frequency w gives a smaller cut- 
off length scale ~ uj~ x I z and leads to an overall shift 
toward the direction opposite to that of the flow (in this 
case, toward the unstable fixed point), and a more broad- 
ened width of the plateau-to-plateau transition. 



V. FLOW DIAGRAMS FOR LARGER a 

So far we have discussed the behavior for a small- 
frequency region, where the dynamical length scale L u ~ 
uj- 1 / 2 behaves an infrared cutoff for the critical phenom- 
ena at the Anderson transition and the dynamical scal- 
ing arguments also hold quite well. It is worth while 
to ask how the behavior of the larger-frequency region, 
typically for uj up to Q.1uj c . looks like in a similar two pa- 
rameter <j xx — a xy plot. Naturally, in the large-frequency 
region, ac conductivities are expected to be no longer 
dominated by the criticality and show a qualitatively dif- 



ferent, rather classical behavior. In this region we should 
adopt raw uj instead of ujL z , because the system should 
be out of the critical (i.e., dynamical-scaling) region. In 
FigHJa), we look at the large-w behavior of the longi- 
tudinal <t xx (u>,Ef = 0) against the frequency oj, from 
oj = O.OIluc to uj = Q.1uj c . In a large-frequency, ac longi- 
tudinal conductivity at the center of LL shows a mono- 
tonic decrease of a xx (e = 0) with uj as consistent with 
Refi22. This clearly signals a deviation from the dynam- 
ical scaling ansatz, which assumes the cjL z -dependent 
conductivity as in eqnJ3] 

The a xx {ui) — <r xy {uj) diagram for the large-w region 
where the flows are indicated with varying frequency 
uj for various values of cf is shown in Fig HJb). The 
w-driven flows show a different pattern from the flows 
in a small-omega region in the previous section, al- 
though it is attracted with decreasing uj to the static- 
QHE fixed points (<J xx ,(J xy ) = (0, ±2)e 2 //i, just as in the 
temperature-driven flows^. This behavior is understood 
that in the large-frequency region, where frequency uj be- 
come comparable to the cyclotron frequency oj c , the fre- 
quency puts a small length scale comparable to the mag- 
netic length, which should naturally induce deviations 
from the critical region and show a pattern of the two 
conductivities different from that in the small-frequency 
region. The existence itself of flows around lu/uj c ~ 1 im- 
plies that the system is not fully dominated by a Drude- 
like behavior, which is consistent with the observation 
of the ac plateau structure in this frequency region, i.e., 
uj ~ 0.1 that corresponds to the THz region. 



VI. SUMMARY 

We have numerically obtained a xx (uj)—a xy (uj) diagram 
for graphene QHE system, and have examined the flow 
for two regimes. In a small-w regime, flows are governed 
by the dynamical length scale posed by the frequency as 
a relevant cutoff length scale for the criticality around the 
Anderson transition. We also discussed a metallic behav- 
ior around the unstable fixed point reflecting the delocal- 
ized state at the LL center. The larger-w regime exhibits 
rather classical flows driven by the bare frequency due 
to the small dynamical length scale comparable to the 
magnetic length scale . 
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